SAYAM Vol-l, Issue-I (2023), Page No-09-19 


as 


—_ ; iy 
wien OD 
epee Mm 


“BES 
fre ss 


Wy} 


ayam 


Journal Homepage: https://sayamjournal.com/ 


Article 


A Non-Linear Type Equation of State and Cosmic Fluid Dynamics 
Shouvik Sadhukhan', Alokananda Kar?, and Surajit Chattopadhyay: 


1Department of Physics; Indian Institute of Technology, Kharagpur 721302, West Bengal, India 


Email: shouvikphysics1996@gmail.com 


2Department of Physics; University of Calcutta; 92 APC Road, Kolkata 700009, West Bengal, India 


Email: alokanandakar@gmail.com 


3Department of Mathematics; Amity University; Newtown; Kolkata; West Bengal, India. 


Email: surajitchatto@outlook.com (corresponding author) 


ARTICLE INFO ABSTRACT 


Keywords: 


Non-linear Equation of 
State, Chaplygin Gas, Van- 
Der-Waals fluid, FRW 
equations, Energy 
Conservation Principle 


Received : 07/09/2022 


Accepted : 02/02/2023 
Date of Publication: 01/06/2023 


In this chapter we have introduced a special type of non-linear equation of state to discuss 
the cosmological evolution mechanism. The new equation of state is a four parameters model 


which can be represented as p = Ap + Bp? — a where B = Af — y. The evolution of universe 


have been interpreted by fluid dynamics. The reconstruction of Chaplygin gas and Van-Der- 
Waals (VDW) fluid equation of states have been done from the parametric analysis of this 
new non-linear model. Different cosmological phases like Quintom, Quintessence and warm 
universe have been discussed here. Finally, we have provided a comparative studies of this 
model with other non-linear fluid solutions. 


1. Introduction 

Recent advancement of experimental cosmology provide 
the proof behind the accelerated expansion of universe. 
Type Ia Supernovae and the CMBR radiation analysis re- 
established the idea of cosmic inflation and late time 
acceleration scheme (Jamil.et.al (2009), Biswas.et.al (2013)). 
Those experimental results brought several questions of 
well known FRW models. The Horizon problem, Magnetic 
monopole problem and Flatness problems are few 
examples of them (Linde.et.al (1982)). For the resolution of 
these problems the cosmic inflation and the late time 
acceleration models have been brought in physics 
(Liddle.et.al (1999), Lucchin.et.al (1985)). The gravity 
provides attraction in local scale. Hence, to bring the 
repulsive nature, the idea of dark energy has got much 
attention. Dark energy is considered to act as exotic cosmic 
system that exerts repulsion against the attractive gravity 
(Copeland.et.al (2006), Steinhardt.et.al (2003), Hughes.et.al 
(2019), Kar.et.al (2021, 2022)). 


There are several dark energy models available in literature 
(Hoyel.et.al (1964), Zeldovich.et.al (1968), Capozziello.et.al 
(2006), Elizalde.et.al (2005)). The most acceptable models are 
Chaplygin gas type fluid dark energy models, scalar field 
type dark energy models and holographic dark energy 
models. There are several varieties of fluid type dark energy 
models viz. the chaplygin gas model (Kotambkar.et.al 
(2021)), modified chaplygin gas (Mamon.et.al (2022), 
Jianbo.et.al (2009)), cosmic chaplygin gas, modified cosmic 
chaplygin gas (Benaoum.et.al (2012), Khadekar.et.al (2019)), 
polytropic type fluid, modified polytropic or Van-Der- 
Waals fluid (Faraoni.et.al (2003), Capozziello.et.al (2002, 
2003), Vardiashviliet.al (2017), Ivanov.et.al (2019), 
Obukhov.et.al (2018), Jantschet.al (2016)) and 
Inhomogeneous fluids (Jamil.et.al (2008), Chakraborty.et.al 
(2008)). These fluid based dark energy models can show the 
self interacting dark matters which can substitute the dark 
energies through Cardassian formalism (Gondolo.et.al 
(2003), Frees.et.al (2002)). The inhomogeneous models can 
establish the theory of dissipation of self interaction of dark 
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matters and can give a pause in the inflationary universe to 


get late time accelerated expansion phase (Kar.et.al (2021, 
2022)). 


The non-linear cosmological models are based on the non- 
linear pressure density relations or non linear equation of 
states. The non-linearity of equation of states are generally 
found in the introduction of self interaction (Nojiri.et.al 
(2004, 2006), Chakraborty.et.al (2019)). The evolutionary 
mechanism which follows those self interactions, are 
generally called Cardassian expansions. Those self 
interactions in Cardassian expansion scheme provides the 
proof of existence of fifth force or fifth dimensions 
(Nojiri-et.al (2004, 2006), Chakraborty.et.al (2019)). The 
decay of those fifth force can bring the graceful exit of 
inflation and start reheating (Albrecht.et.al (1982)) for late 
time accelerations. The non-linearities in equation of state 
are basically a combinations of linear dark matters and 
some non-linear fifth force based potentials which are 
active in Hubble scale but negligible in local scale. This 
force causes accelerated expansion against the 
gravitational collapsing (Kar-.et.al (2021, 2022)). 


In this present work we have used a new type of four 
parameters equation of states. This equation of state can be 
a generalization of both Chaplygin gas (Jianbo.et.al (2009)) 
and VDW fluids. This fluid model can be derived using 
binary type self interaction between linear dark energy, 
polytropic fluids and chaplygin gas (Elizalde.et.al (2018), 
Brevik.et.al (2018), Wuet.al (2008)). We have only 
discussed the fluid dynamical properties and its effect on 
universe expansion into this paper. We have derived 
evolution scheme of the energy density and pressure of this 
model under Einstein gravity with non changing 
gravitational and cosmological constants. Finally, we have 
re-established the chaplygin gas (Benaoum.et.al (2012), 
Pourhassa.et.al (2014), Kotambka.et.al (2021)) and VDW 
fluid from this model along with the discussion of some 
cosmological phases. 


The paper is oriented as follows. In section 2 we have 
discussed the mathematical basis of our work. In sections 3 
and 4 we have given the reconstructions of different fluid 
models and different cosmological phases respectively. 
Finally, in sections 5 and 6 we have provided the results 
analysis and concluding remarks of our work. 


2. Mathematics on the New Non-linear fluid 


In this present work we want to study a new type of non- 
linear fluid model which can replace linear dark matter to 
discuss the expansion scheme of universe. We have used 
the action of constant gravitational constant with absence 
of cosmological constant type dark energies. We have 
taken 87G = c = 1. The necessary action for our model can 
be written as follows. 
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oS | d*x(,/—gR + Ly) 


Now this action has been used on the isotropic and 
homogeneous geometry i.e., the FRW geometry. The line 
element has been taken as follows. 


ds* = dt? — a*(t)(dr? + r7dQ?) 


The line element represents a flat universe and dQ? = d@? + 
sin* 6dp*. Now using least action principle on the action 
given above, we can find the following dynamical equation 
of cosmic mechanics i.e. the Einstein field equation. 


KT, 


1 
Ruy = 5 Iuvk = Klyy 


Here, k = 1 in SI unit. The term T,,, comes from the cosmic 
matter Lagrangian L;. The relation can be explained as 
follows. 


=, ttl aos) 
fag og! 
The tensor representation of the energy momentum tensor 


can be written as follows using the knowledge of isotropic 
energy distribution. 


pv 


Thy = (py + PF )uytly —~ Pf Inv 


Now, using the isotropic and homogeneous line element 
into the Einstein field tensor, we can find a couple of 
dynamical equations that can discuss the evolution of 
isotropic and homogeneous universe. The couple equations 
are called FRW equations which are as follows. 


3H? = Pf 


3H? + 2H = —py 


For any cosmological models, the fluid energy must be the 
conserved, i.e energy momentum tensor must be divergence 
less. Similarly, the geometry must follow the divergence less 
condition also. Hence we must have to satisfy V,G,, = 
VuTyy. Thus, we found the following energy momentum 
conservation equation. 


py + 3H (pf + Pr) =0 


Now the solution of cosmic matter energy density and 
pressure should follow the energy conservation equation. 
The energy momentum conservation equation can provide 
the scale factor dependent variations of energy density and 
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pressure which again can provide the time variation of 
scale factor by substitution of energy density and pressure 
from FRW equations. During the derivation of scale factor 
variations of energy density and pressure we must discuss 
the equation of state or pressure density relations of our 
model. This is given as follows: 


Pam 
Dp = Ape + Boz ——~| 
PF 


The non-linearity of the equation of state can provide 
dimensionally non unit parameters. Here, the dimensions 
of A,B and C are [1],[p]~* and [p]%*? respectively. The 
parameter B can be dissociated as B = AB — y where the 
dimensions of 6 and y are same as B. Now using the 
equation of state into the energy momentum conservation 
equation, we can find the following integral. 


a 
d 
| ieee +C,=-—3lna 


Bpet? + (A+ 1)p¢** —-C 


We have to find the solution of this equation which 
depends on four fluid parameters A, 8B, y,C and one power 
parameter a. Hence, we have applied the fluid dissociation 
scheme to solve the energy density and pressure of our 
model. In the fluid dissociation method we have 
dissociated the fluid equation of state into two parts. One 
is Dr, = Apy, — a and another one is pg, = (AB — Y)PF,- 
1 
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Here, the first part of equation of state act as cosmic fluid on 
p-brane whereas the second part of EOS act as the D-brane 
or fifth force based self interaction between the cosmic fluid 
particles. This dissociation is viable in physics as this second 
component works as fifth force to modify the effective 
pressure but don't work as second fluid system. Hence, the 
energy momentum tensor can also be dissociated into two 
equations without having any interaction scenario. The self 
interaction can not be shown with usual interaction scenario 
available in cosmology literatures. Hence, the effective 
energy density and pressures can be written as follows. 


Pr = Pr, + Pp 


C 
Dp = Pr, + Pp, = App, + (AB — ) pF, — oF 


fi 


The energy momentum conservation equations for the 
dissociated system can be written as follows. 


pp, + 3H (py, + Pp) = 0 


Pp, + 3H(pp, + Pp) = 0 


Now we solve those energy densities and pressures as a function of scale factor with usual ways. Hence, The solutions 


are as follows. 


1 


C CG \re@ 
Pf, = (; +A a aria 
_ 1 
Pf ~ Ca — (AB—Y) 
1 a 
G G. Tra C G “Ta 
Df ==A (; +A i “aa —C¢ (; +A + ae 


Pp = (AB -y)( 


1 
C,a* — (AB — 5) 


2 


Hence, the effective energy density and pressures can be written using the above relations. The effective density is as 


follows. 


c @ 
Crs (; +A + qadtAy(ata) 


1 
ae 1 


* Ca? — (AB —1) 
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And, the effective pressure can be written as follows. 


1 1 
- Cc Cy The Cc Cy ra 1 - 
se (& * Dy, seat ees) AB SY) (aa) 


Hence, the effective equation of state parameter can be written as follows. 


1 1 
G. Cy T+a C Cy “Ta 1 2 
ae (( tat poraa) -Crgat pada) "+ 48-1) (<2 =p = >) 


1 


PF G 1 


( C 4 ae 4 
1+A g3QtA4)(1t+a) C,a3 = (AB = Y) 
Now using the above equations and controlling the parameters A, 8, y,C and a we can reconstruct the chaplygin model 
solutions as well as VDW fluid solutions for cosmic evolution. Different phases and their smooth transformations through 
entire cosmic era can be discussed using the values of equation of state parameters. 

3. Generation of different fluids 
This new model discussed in the above section can be used to reconstruct different cosmological non-linear models with 


some parameter analysis. We can discuss the Chaplygin gas, Modified Chaplygin gas, Van-Der-Waals fluid and 
Polytropic fluid. The models can be reconstructed as follows. 


3.1. Chaplygin Gas 


For reconstructing the equation of state of Chaplygin gas, we have to choose the parameters as A — 0, (AB — y) > 0 and 

a > 1. Hence, the equation of state should be as py; = — <, pressure, density and EOS parameter wy can be written as 
f 

follows. 


Cy “2 
a (| 


and 


oy 
Cy 
( qeo(1+A4) ) 
Hence, these results can produce pure Chaplygin gas and different cosmic phases from it. 


3.2. Modified Chaplygin Gas 


Now we want to reconstruct the equation of state of Modified Chaplygin gas from our model. For this purpose we must 
have to consider (AB — y) > 0. Hence, the equation of state can be taken as pr = Apr — a The pressure, density and 
f 


EOS parameter w, can be written as follows. 
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1 


4 & Cy ita 
Pf = (; LA + aa] 


1 a 
C C; ita € C, “ita 
de (; A ars] so (; cA iar 
and 

C C lS C C a 

Co en ee 1+a oe on no 1+a@ 

A (( oa aura) C(qpat aE) 
wm 7: 000... Seer 7 
oc <2 Gy, e- <2 1+a 
(7 +At aura) 


3.3. Van-Der-Waals Fluid with |Bp;| < 1 


The reconstruction of Van-Der-Waals fluid can be done using the assumptions of parameters as C > 0. Hence, the 
equation of state should become pr = Apr + (AB — y) p; which is a special form of VDW model with condition of 

|Bp,| < 1 (details is given in appendix). The mathematical solutions of energy density and pressure as well as equation 
of state parameter can be written as follows. 


1+A 
Pr (= AB) + pro + MGA 
1+A 1+A m 
a & AR) A papi nea) Tap) —AB) + Pmo(l + ee) 


and 


1+A 1+A ; 
4(Gaapeacdeaars) + 6-9 (Gaara) 
a a 
(y — AB) + pro (1 + A)a3@*4) 


This model can discuss the cosmic phases and their evolutions with three parameters. For this model we found some 
constraints on the parameters and scalar factors which can be explained as follows. We can assume that the scale factor 
can run only within the limit (see Appendix) 


(f tap == ae Pe (S PA or y= ee 
Pro(1 + A) Pro(1 + A) 


Similarly we can provide the limits of pressure of the activity of that equation of state as follows. 


4, AB), 4 B=”) 
Bp BP 


Hence, pressure can be negative with the dupplet as (vy, f) € (vy << 0,8 >>0)U(v >> 0,6 << 0). 


3.4. Polytropic Fluid 


For the reconstruction of polytropic fluids using our present model, we have to bring the constraints as A > 0 and C > 0. 
Hence, the equation of state should become py, = (AB — Y)p?- The solutions can be written as follows. 
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1 
*r C,a3 — (AB -7) 
1 2 
= 8 - (Car aE=D) 


and, 


Wr = 


1 


C,a° — (AB —y) 


These solutions can provide an universe model which follows quadratic fluid formalism. 


4. Different Cosmological Phases and their fluid 
dynamics 


We aim to discuss different cosmological phases with 
parametric analysis from the new non linear model of our 
work. The model and its solutions are dependent upon four 
equation of state parameters and one power parameter. 
They are A, B,y,C and a respectively. Now the magnitudes 
of those parameters can provide different cosmological 
phases, their evolutions and smooth phase transitions. 
Here, we have discussed the parametric analysis for getting 


some cosmic phases and their transformations in table 1 . 
The table 1 represents the region of the magnitudes of the 
parameters which can be used to discuss different cosmic 
phases and their transformations as well as their evolutions. 
We have shown the densities and pressures for those 
phases graphically in next sections. We have chosen the 
values of the parameters to plot the figures according to the 
region derived in table 1. a has been considered positive to 
bring the end of each phase with evolution w.r.t. scale 
factor. 


Phantom Phase 


Quintom Transformation 


Quintessence Phase 


CDM Dominant Phase 


WDM Dominant Phase with Radiation Region 


A<0,B<OandC>0 
hence, B > 0 and y > 0 witha <0 


A>0,B>OandC >>0 


hence, either B >> 0 ory << 0 witha >0Oora<1 


A>0,B>OandC >>0 
hence, both 6 >> 0 and y << 0 witha >> 0 


A=0,B=OandC =0 
hence, AB = y with a = arbitrary 


A>0,B>OandC << 0aswellasC <0 
as well as C < 0 
hence, both 6 >> 0 and y << 0; or also 
hence, either B >> 0 ory << 0 witha > 0 


Table 1: Parametric Analysis of different Cosmological Phases 


5. Result Analysis and Comparative studies 


We have used pre predicted magnitudes of the parameters 
and our results is capable of establishing the Phantom 
phase, Quintom region, Quintessence phase, Cold Dark 
Matter or CDM phase as well as Warm Dark matter i.e. 
WDM and radiation properly with exact values of the 
equation of state parameters. The results has been given in 
figure 1 from a to o. We have given the results for different 
phases in different graphs. All the densities are positives in 


14 


magnitudes and pressures are negative except for the 
radiation phase. For Phantom we have found w, < —1, for 
quintom its around w, ~ —1. For Quintessence and CDM 
the equation of state parameter gave 0 > we > —1 and wy = 
0 respectively. For radiation, both pressure and density 
provides positive values and hence, the equation of state 
parameter comes as wy ~ 0.33. All these results found are 
in their evolving state i.e. all the energy densities and 
pressures are changing with scale factor. Hence, we could 
provide smooth transitions during the phase transitions of 
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universe. In our analysis, we consider only the fluid 
dynamics and our results are in agreement with the 
experimental results. Except the Phantom dominant 
universe, we found lpr| > pl. Hence, these models are 
stable with Null energy conditions, Weak energy 
conditions as well as Dominant energy conditions in the 
perspective of fluid dynamics. All of these models provide 
|p¢| << |3p,| and thus, with the violation of strong energy 
conditions they can provide accelerated expanding 
universe solutions. The results shown in figure 1 provide 
specific region of the scale factors where they can only 


Cosmic Energy Density vs Redshift 


Cosmic Pressure vs Redshift 


Sadhukhan, S., Kar, A., and Chattopadhyay, S..,2023 

provide the cosmological phases discussed here. These 
phenomena strictly proved that no cosmic phases discussed 
here, can be stable for permanent period of time. All of them 
must start and end through some phase transitions. These 
phase transitions should be controlled by the four cosmic 
parameters and one power parameters discussed in our 
model. Thus, we can conclude that the parameters of our 
model are also time varying during the phase transitions 
but constant or almost constant during specific phases. The 
change of the magnitudes of those parameters provide an 
idea about the order of phase transitions for those phases. 
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Figure 1: Black for Radiation and WDM, Blue for CDM, Red for Quintessence, Pink for Quintom and Purple for 
Phantom universe. 
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Cosmic Energy Density vs Redshift 
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Figure 2: Black for Radiation and WDM, Blue for CDM, Red for Quintessence universe 
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6. Conclusion 


We have finally, established the new non _ linear 
cosmological fluid model through the detailed discussion 
of fluid dynamics. The whole discussion proceeds with the 
scale factor variation of energy density and pressure for a 
new non linear equation of state under Einsteinian gravity 
action with non varying gravitational constant. We have 
matched our results with the well established values of the 


Sound Speed Velocity vs Redshift 
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equation of state parameter for different cosmological 
phases. Our results also violate strong energy condition to 
provide acceleration in the expansion scheme of our model. 
The null energy condition for all those phases is satisfied, 
because we found |p,| > |p,|. Hence, phenomenologically, 
our model obeys the first law of thermodynamics beyond 
phantomic region of universe. The details of the 
thermodynamical analysis is beyond the scope of this work. 


Deceleration vs Redshift 


az) 


(b) q(z) vs z 


Deceleration vs Redshift _ 


az) 


-21.0 


-21.2 


(d) g(z) vs z 


Figure 3: Pink for Quintom and Purple for Phantom universe 


The variation of deceleration parameters and sound wave 
velocities have been presented in 10 graphs of figures 2 and 
3. Those graphs have been given to represent the stability 
of our model discussed in this paper. The sound wave 
velocity should be positive when the model is stable under 
scalar perturbation. Hence this four parameter equation of 
state based cosmic fluid model is stable under perturbation. 
The deceleration parameters for each phase provide 
negative values which proved that all those models are 
representing cosmic acceleration. 
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Appendix 


The generalized Van-Der-Waals model or VDW fluid 
equation of state for cosmic dynamics can be written as 
follows. 
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Now this equation of state can be expanded into binomial series as follows. The denominator of the first term of RHS of 


the above equation has been expanded with binomial series. 


Pr = Apr (1+ Bpr + (Bp,) + (Bos) +) - 107 


Now, if we consider |B p,| < 1, we can converge the expanded equation of state as 


pr = Ape + (AB —y) pF 


Now in presence of that condition, we can provide the following analysis which can control the cosmic dynamics under 
that equation of state. The condition |B p,| < 1 can be extended into following inequality. 


Bp;|<1>-1< Bp, <1 

Hence, we can provide two concluding statements which are as follows. The mathematical form of the first statement is 
as follows. 

1 Sneed 1 

—_— p a 

po’ “B 
and, the second statement is as follows. 

ie 2 1 

Pf Pr 


The first statement (equation 40) represents the boundary 
condition of energy density. For having fixed parameter f 
there must be some upper and lower boundary for energy 
density in the cosmic fluid system under the specific type 
of non-linear model or specific VDW model. Beyond this 
limits the cosmic fluid will get some phase transition and 
hence it starts to follow different equation of states. 


The second statement provides the nature of variation of 
parameter £. According to the representation the parameter 
f should vary hyperbolically with energy density with a 


limit of (-4.4), Whenever the parameter crosses the 
f Pf 
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